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Abstract 

I have performed experimental measurements of 
acoustic emission signals resulting from the drying 
process of Phyllostachys Pubescens bamboo. The em¬ 
phasis was on identifying individual events, and char¬ 
acterize them according to their time span and energy 
release. My results show a histogram of experimental 
squared voltage distributions nicely fit into a power 
law with exponent of —1.16, reminiscent of scale free 
phenomena. I have also calculated the average sig¬ 
nal shape, for different time spans of the system, and 
found an asymmetrical form. The experimental evi¬ 
dence points to the system having an isolated large 
crack at the beginning of the simulation. 

1 Introduction 

Bamboo is a giant grass, with remarkable mechani¬ 
cal characteristics Q Nonetheless, as it is a natural 
material, its drying requires special care in order to 
prevent the cracking due to dessication [5,. In some 
aspects Bamboo drying is similar to wood drying: 
there is a first fast drying stage, in which free water 
is removed from the plant; followed by a slow stage in 
which bond water, from within the cell walls, has to 
be removed m- The main difference between drying 

1 Values of the elastic modulus of the fiber cell wall of ei¬ 
ther 10.4 GPa [13] or 16.1 GPa [12|have been reported in the 
literature 


this two natural materials comes from the structure 
of their tissues. Unlike wood, bamboo has a special 
direction, the culm, along which different tissues are 
aligned: its fibers, conductive material and parenchy¬ 
matous tissues 0. 

Very few research groups have studied the relation¬ 
ship between shrinkage and fracturing in bamboos. 
Overall shrinkage of bamboo canes can be attributed 
to the shrinkage of the fiber-phloem bundles, as well 
as to the shrinkage of parenchyma tissue. In the lat¬ 
ter, the collapse of cells can be an important factor. 
Obataya et. al. have shown how the slow drying pro¬ 
cess can lead to the collapse of Arundo donax cane, 
HUS]. Those authors claim that an slow drying sched¬ 
ule increases the intensity of collapse. 

On a more general fracture case, Habibi and Lu 
have studied the path followed by different cracks 
within the bamboo culm, by means of micro inden¬ 
tation experiments as well as tensile tests [2]. In¬ 
terestingly, the micro crack grows mainly in the in¬ 
terface between the parenchyma cells. When the 
crack reaches a fiber it is sometimes possible for it 
to continue within the fiber, again in the interfaces 
between fiber cells. They have also shown how the 
voids within the bamboo have an effect of crack grow 
deflection and crack tip energy absorption. 

I would expect a different path of growth for 
dessication induced fractures compared to indenta¬ 
tion; as drying stresses are more homogeneous, not 
having a unique defined direction of deformation. 
My attempts to look for the fracture in the surface 
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proved futile; which may be caused by fractures being 
smaller than the resolution of my optical microscope, 
or by them being inside of the samples themselves 
[5|. Then, the problem calls for a different way to 
characterize the fracture evolution during drying, so 

1 turned the research focus into acoustic emission. 
As fractures are created at every scale, some of the 
energy is released in the form of an acoustic wave. 

In the present letter I have characterized the dessi- 
cation induced fractures in bamboo Phyllostachys 
Pubescens. To do so, I performed drying experiments 
in which I recorded the acoustic signal produced by 
the sample. Since it was not subject to any other 
strain, it is safe to assume that the sounds are gener¬ 
ated by the dessication induced strains. The goal of 
this measurement is to obtain information about the 
drying induced fractures. The hypothesis being that 
the average shapes of fractures, the distributions of 
energy released in an avalanche event, and also the 
waiting times between avalanches; all have useful in¬ 
formation about the process of the drying induced 
fracturing. 

The organization of the paper is as follows. In 
section 2 I describe the experiments, materials and 
methods used. In Section 3 I will present my results 
and conclusions. 

2 Materials and Methods 

Two kinds of experiments were performed, namely, 
on characterizing the process of drying and on char¬ 
acterizing the acoustic emissions. For the former I 
measured the weight as a function of drying time. 
For the latter, I measured the acoustic emission re¬ 
sulting from the drying process. 


Label 

Number 

Samples 

of 

Length 

Sample 0 

4 


1 cm 

Sample 1 

1 


10 cm 

Sample 2 

10 


10 cm 


Table 1: Phyllostachys Pubescens inter-node samples 
used in the present study. 


2.1 Materials 

Bamboo samples were prepared from a 4 year old 
culm Phyllostachys Pubescens , harvested on June 
2011, from Friedrich Eberts’ plantation in Chiavari, 
Italy. Three sets of inter-node cut samples were la¬ 
beled and prepared, as shown in the table. On Sam¬ 
ple 0, I performed drying experiments aimed to char¬ 
acterize the likelihood of collapse and fracture. On 
Sample 1, I took the weight measurements during 
drying, in order to characterize the speed of dry¬ 
ing. Finally, on Sample 2, I measured the acoustic 
emission of breaking events during drying. Between 
measurements the samples were kept on a climate 
controlled room, with relative humidity of 80%, to 
minimize the amount of drying happening before the 
oven treatment. I performed the drying experiments 
using a SalvisLAB Pantatherm D oven. 

2.2 Methods 

Two sets of experiments were performed. First, for 
the characterization of the speed of drying, 20 inter- 
nodal bamboo samples of 10 cm width were dried in 
the oven at a set temperature of 100 C. In this ex¬ 
periment the percentage of weight loss of the bamboo 
samples as function of time was monitored; and later 
averaged to build the plot. 

The second set of experiments characterized the in¬ 
tensity of the microscopic cracks by measuring their 
acoustic emission during drying. The drying temper¬ 
ature was again 100 C, and the digitization was per¬ 
formed using a standard condenser microphone and 
a digital sound card at a rate of 9600 bps. The volt¬ 
age on the microphone is proportional to the energy 
released as acoustic waves. After the recording, the 
whole data set was screened to check for noise coming 
not from the drying specimen. 

The recording took place during the first 3 hours 
of drying. Background noise, both acoustic and from 
the electronics, was filtered by means of setting a min¬ 
imum threshold for recording a voltage signal (See 
Fig. 0 - Experiments were performed overnight, to 
reduce the amount of noise from the environment. 
Later on, the recordings were split in 6 s samples; I 
listened to each one of them in order to exclude the 
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sections of the recording that presented strong noises, 
probably arising from the oven. Finally, the signals 
were converted in text files. 

Now, I claim that the continuous acoustic emis¬ 
sion events can be interpreted as produced by the 
breaking of bamboo structures; be them fibers, con¬ 
ductive materials or parenchyma. Therefore they are 
be called indistinctly crack events or acoustic events. 
A single event is defined as a continuous signal that 
lies above the minimum recording threshold. 

Since my objective is the characterization of the 
fracture process, the exact relationship between the 
voltage measured by the microphone and the ampli¬ 
tude of the original sound wave in dB does not need 
to be explicitly known. Moreover, this calibration is 
a quite complex process, that involves the electronics 
of the sound card and the software of the computer 
used to the recording (OS X, a proprietary and closed 
source operating system by Apple Inc.). Therefore I 
will report my experiments in arbitrary units of am¬ 
plitude: a.a., and its square is in arbitrary units of 
energy, a.e.. 


3 Results and Discussion 

3.1 Speed of drying experiments 

The average of the weight of the bamboo samples as 
function of drying time is shown in Fig. [2] Let me 
focus first on the black dots. It seems like there is a 
linear trend in the semi-log graph for the first three 
hours of drying. This can be related to a rate process, 
which can be the removal of the free water (filling the 
interior of the cells, not chemically bounded). This 
is followed by a slowest drying phase, between hours 
three and eight. This may be due to the removal of 
the bound water (from inside the cell walls, chemi¬ 
cally bounded). The thicker samples (red diamonds) 
only show the first stage, but the moisture content 
has not fallen below 2%, so there is still bond water 
in the sample. 


3.2 Energy release histograms and 
waiting times 

The study of fracture phenomena can be tackled from 
different disciplines. From the point of view of crit¬ 
ical phenomena and statistical mechanics, fracture 
is seen as a series of bursts; whose size and tempo¬ 
ral evolution characterizes the phenomena. Recently 
Papanikolau et. al. have taken the Barkhausen noise 
and studied the functional form of the noise emit¬ 
ted by the avalanches 0. Besides the signal aver¬ 
age shape, they present voltage distributions, power 
spectrums, and distribution of avalanche sizes S and 
durations T. Their avalanche size histogram shows a 
power law distribution at small sizes, followed by a 
tail. 

In the present case, Fig. [3] shows the probability 
distribution of the square of the voltage recorded by 
the microphone, which is proportional to the energy 
liberated in the fracture. A power law distribution 
of avalanche sizes is one indicative of scale free be¬ 
havior, as seen for instance in the Gutenberg-Richter 
law of avalanches, and in the Barkhausen noise; even 
though the mechanisms that generate the fractures at 
different scales, there is not a single size of fracture 
that rules the breaking phenomena. In this case the 
distribution of the square of the sizes follows a power 
law (which is proportional to the energy released in 
the fracture event), 

P(A) oc A~ a , (1) 

with an exponent of a = —1.16. It is interesting 
to notice that the power law seems to fit for three 
decades. In the present system small fractures and 
large ones correspond to different kind of tissues, and 
can be caused by accumulation or by direct strain; so 
it is not at all self evident that the fracture distribu¬ 
tion has to be a power law. 

A second characteristic measurement that can give 
insight into avalanche processes is the waiting times 
between avalanches. For earthquakes this distribu¬ 
tion is known to follow a unified scaling law, mostly 
a power-law with a decaying exponential tail PQ. In 
the present study, on the contrary, waiting times be¬ 
tween avalanches show an exponential decay behav¬ 
ior; as can be seen in Fig. [4] In this case the time 
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constant is 0.0574 s. The difference with the earth¬ 
quake case can be explained from the fact that in this 
system the allowed number of cracks is reduced as 
there are a finite number of elements to crack or de¬ 
laminate; while for earthquakes the possible number 
of rearrangements and fractures is orders of magni¬ 
tude larger. In this way, this system is closer to the 
granular materials. In a recent paper, Michlmayr and 
Or have studied the relation between the grain-scale 
mechanical interactions in sheared granular materials 
and the generated acoustic emission characteristics 
[3] . Their implementation of a damage accumulation 
fiber bundle model (similar to the one by Kun et. al), 
shows an exponential decay of the number of events 
as function of the energy (their Fig. 3, inset); which 
is the analogous of the waiting times for that system. 

3.3 structure of energy release 

To further characterize the structure of the 
avalanches, I want to know when do the big 
avalanches occur with respect to the start of the dry¬ 
ing process; as well as whether large energy release 
indeed means long avalanches. For the former ques¬ 
tion, to identify the evolution of energy release, the 
size of the avalanche is averaged over 5 second sound 
intervals and plotted as function of the number of the 
sound interval; and can be seen in Fig. [5] It is pos¬ 
sible to identify two stages, early and late, separated 
by the sound interval number 1000. In the first stage 
the intensity is a large and decreasing quantity. At 
the late stage, there is a steady increase of the size of 
the avalanches as a function of time. The black lines 
represent the error bars. 

A possible explanation for the shape of this graph: 
large values for short times, followed by a steady in¬ 
crease is this: Initially, the drying process generates 
large fractures; creating free surfaces in a fast pro¬ 
cess. The presence of this free surfaces hinders the 
appearance of large cracks. After this follows a sec¬ 
ond phase, in which further fractures may happen 
mainly for stress accumulation. This two behaviors 
are seen in other computational models nxa • In that 
previous work, that I did with other collaborators, we 
shown how the statistics of bursts of breaking events 
in the case of shrinkage of a thin layer of 2D material 


could suggest the presence of an underlying critical 
point. In that case, the two states are related to the 
way the fracture is created, either by coalescence of 
micro-cracks or by the appearance of a large isolated 
crack that spanned the sample. 

To show the relationship between the size of the 
avalanches and its duration, I present the correla¬ 
tion between the amount of released energy in an 
event and its time span in Fig. [6] as an scatter plot. 
The lower cutoff comes from the imposed threshold of 
measurement. The discreetness in time span values 
is a consequence of the sound sampling rate of the 
acquisition transducer. Clearly, there is more disper¬ 
sion on the crack events that take longer and release 
a larger amount of energy, while the distribution nar¬ 
rows for shorter times. Furthermore it is clear that 
the relationship between energy release and time span 
can be fitted to a power law spanning several decades. 

3.4 Shape of the acoustic emission 
events 

A final way to characterize the process is to study the 
average shape of the acoustic emission signal, for dif¬ 
ferent sizes. Knowing which shape do small and large 
avalanches have can be useful to identify the kind of 
process that generates this fractures. As a matter 
of fact, in a previous paper [5] showed a temporal 
average avalanche shape, rescaled to unit height and 
duration, has a parabolic shape for small avalanches, 
that flattens towards a unit pulse shape for longer 
avalanches; in the case of Barkhausen noise. 

To define the intervals that correspond to 
avalanches of similar size, I performed an histogram 
of the duration of avalanche events, Fig. [7] The 
number of avalanches varies over 4 orders of mag¬ 
nitude. In view of this, I hand picked six inter¬ 
vals: At e (0,0.002)s,At e (0.002,0.005)s,At e 
(0.005,0.01)s,At e (0.01,0.012)s, At e (0.012,0.02)s 
and At £ (0.02, 0.04)s. Roughly speaking the first 
three correspond to the most populous bins (10 2 to 
10 4 fractures), short avalanches; and the latter to the 
large avalanches, more than 0.01s. 

Fig. 0 shows the average avalanche shape for the 
different bins, as a function of the normalized time 
span. In the main plot I show the first five averaged 
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signals, with the longest avalanches are shown, but 
faintly; while in the inset all the six averaged signals 
can be seen. Unlike the Barkhausen effect case from 
the literature, 13J in which there is clear symmetry; 
in this case most of the signals are not symmetrical: 
there is a large peak followed by a decay. This can be 
seen in all intervals but one (At £ (0.012, 0.02)s). In 
my view, the asymmetrical shape shape hints for the 
phenomenon to be that of a single breaking event, 
and the decay in the signal corresponds to attenu¬ 
ation of the noise. Short signals would correspond 
to small (or weaker) elements, as parenchymatous 
tissue; while long signals to larger (or stronger) el¬ 
ements, as fiber cells. Going back to the average 
signal in the range At £ (0.012, 0.02)s, the symmet¬ 
rical shape here would pinpoint to a different mech¬ 
anism, namely the correlated breaking of similarly 
sized elements; a self induced growth process in which 
a breaking triggers similar events. In terms of the en¬ 
ergy, in Fig. m it can be seen that the average shape 
of the largest avalanches has a peak of energy dissipa¬ 
tion which is almost two orders of magnitude larger 
than that of the medium and small size avalanches. 
Would be interesting to test whether this is caused 
by the composite nature of the bamboo, or is rather 
the result of the dynamics of the avalanche forma¬ 
tion; and would occur also on a more homogeneous 
material. 

4 Conclusion 

The experimentally measured acoustic emission sig¬ 
nals resulting from drying of bamboo Phyllostachys 
Pubescens shows evidence of free-scale phenomena. 
Firstly, the histogram of experimental squared volt¬ 
age distributions nicely fit into a power law with 
exponent of —1.16; while the experimental waiting 
times do not follow a power law, rather an exponen¬ 
tial of slope —17.4 possibly due to the finite size of 
the system. Secondly, the scatter plot of energy as 
function of avalanche time span also shows a power 
law that spans over 4 orders of magnitude. I have 
also studied the average avalanche shape, and found 
that both short and large avalanches have an asym¬ 
metrical shape with a peak at the beginning. This 


implies that for those cases there is mainly a single 
crack event. For medium sized fractures, on the other 
hand, I have found a more symmetrical shape; which 
could be attributed to fracture growth by a self rein¬ 
forced process. 

The main question that rests to be answered is 
whether this behavior is particular to the drying of 
bamboos; or can be found in other materials. This 
can be studied easily from numerical simulations; 
which is what I plan to accomplish in the future. 
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Figure 1: (Color Online). Examples of voltage on the 
microphone, blue continuous line (proportional to the 
sound intensity), and set noise threshold. (Upper) A 
cascade of small sized breaking. (Lower) A single 
large crack. (Sound included as supplementary ma¬ 
terial.) 



Figure 2: Moisture content as a function of temper¬ 
ature for: one 10 cm sample (diamonds) and four 1 
cm samples (black point at the average value). 



Figure 3: Experimental squared voltage size distri¬ 
butions N( A). The line drawn has a slope of —1.16. 
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Figure 4: Experimental waiting times between con¬ 
secutive avalanches. 
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Figure 6: Experimental scatter plot of energy as func¬ 
tion of avalanche time span. 
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Figure 5: (Color online) Experimental averaged en¬ 
ergy release per avalanche as function of time. I only 
show the upper error bars. 
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Figure 7: Histogram of duration of avalanche events 
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Figure 8: Normalized averaged signal shape, for 9 dif¬ 
ferent bins of time lengths. 



Figure 9: Normalized average square signal shape,for 9 
different bins of time lengths. 
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